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We report the design and characterization of a self-assembled, locally resonant acoustic metama¬ 
terial for Lamb waves, composed of a monolayer of 1.02 pm polystyrene microspheres adhered to a 
1.3 pm thick free-standing silicon membrane. A laser-induced transient grating technique is used 
to generate Lamb waves in the metamaterial and measure its acoustic response. The measurements 
reveal a microsphere contact resonance and the lowest frequency spheroidal microsphere resonance. 
The measured dispersion curves show hybridization of flexural Lamb waves with the microsphere 
contact resonance. We compare the measured dispersion with an analytical model using the con¬ 
tact resonance frequency as a single fitting parameter, and find that it well describes the observed 
hybridization. Results from this study can lead to an improved understanding of microscale contact 
mechanics and to the design of new types of acoustic metamaterials. 


Locally resonant acoustic metamaterials are a type of 
composite material that have been the subject of intense 
study over the past fifteen years 1] due to their ability to 
exhibit extreme [2], anisotropic j3], negative ![4|, strongly 
absorbing @], and locally-tailorable [q] effective proper¬ 
ties. These unique properties stem from the interaction of 
propagating acoustic waves with subwavelength resonant 
elements forming the composite jY]. This can present fab¬ 
rication challenges as wavelengths are reduced, includ¬ 
ing challenges involved with fabricating the resonators 
themselves, and in fabricating large areas of composite 
containing such complex microstructure. In this respect, 
colloidal self-assembly is a promising solution, as it has 
been shown to enable the simple, inexpensive, and fast 
fabrication of complex, ordered structures composed of 
nano- or microscale elements in one to three dimensions 
[9]. These advantages have driven the use of self-assembly 
strategies in multiple areas, particularly in the design of 
phononic crystals EHEUi photoni c cry stals EMI , plas- 
monic sensors and nanostructures [13| , and surfaces with 
tailored wettability E3 ■ Despite this wide use, there re¬ 
main few examples of locally resonant acoustic metama¬ 
terials fabricated using self-assembly techniques EE EH ■ 

Lamb waves are a type of acoustic waveguide mode 
that occur in thin clastic plates and membranes E3 that 
play an important role in nanomechanical resonator Ell 
and sensing applications [lfij . They have also been uti¬ 
lized in studies of sub-THz phonon transport with impli¬ 
cations for the understanding of nanoscale thermal phe¬ 
nomena [2(| ■ Recently, several examples of locally reso¬ 
nant metamaterials for Lamb waves have been explored 
in both theoretical [2l| and experimental |22| settings. 
However, in each of the experimentally realized cases, 
the resonant elements forming the metamaterial had fre¬ 
quencies ranging from the audible regime to a few mega¬ 
hertz and were fabricated using conventional machining 
or microfabrication techniques. 

In this work, we present the first realization of a 



FIG. 1. Illustration of the metamaterial self-assembly fabri¬ 
cation procedure, including: (a) monolayer formation at the 
air/water interface; and (b) transfer of the monolayer to the 
membrane, (c) Representative image of the polystyrene mi¬ 
crosphere monolayer. Schematic illustrations of the (d) theo¬ 
retical model and (e) experimental setup. 


self-assembled, locally resonant metamaterial for Lamb 
waves. The metamaterial is composed of a self-assembled 
monolayer of microspheres adhered to the surface of a 
thin silicon membrane. As a result of the adhesion, the 
microspheres have a contact resonance where the micro¬ 
sphere oscillates like a rigid body, with a localized re¬ 
gion of elastic deformation around the point of contact 
with the substrate. In our system, the contact reso¬ 
nance plays the role of the metamaterial locally reso¬ 
nant element and has a frequency of 200 MHz. This 
type of metamaterial may have potential future advan¬ 
tages for acoustic wave tailoring applications, as arrays 
of macroscale spherical particles have been shown to sup¬ 
port unique nonlinear dynamical phenomena due to the 
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Hertzian relationship 23| between spherical elastic par¬ 
ticles in contact [2J, [23 ■ To characterize the metamate¬ 
rial, we use a laser-induced transient grating technique 
[26l l27| to excite long-wavelength (relative to the parti¬ 
cle size) Lamb waves in the metamaterial and measure 
its acoustic response. This approach was recently used to 
study the interaction of microsphere contact resonances 
with Rayleigh surface acoustic waves (SAWs) [16] . In 
the dispersion curves derived from our measurements, we 
identify both the presence of the lowest frequency micro¬ 
sphere spheroidal resonance and hybridization phenom¬ 
ena typical of locally resonant acoustic metamaterials, in 
the form of a previously unobserved phenomenon: the 
hybridization of fundamental flexural (Ao) Lamb wave 
modes with a contact resonance of the microspheres. We 
develop a model that assumes an axial contact resonance 
and yields an analytical expression for the dispersion re¬ 
lation. In the calculated dispersion curves, we observe 
that the microsphere contact resonance not only couples 
with Aq modes, but also with the fundamental dilata- 
tional (So) Lamb wave modes, which results in coupling 
between flexural and dilatational modes. We find good 
agreement between our measurements and the model us¬ 
ing the microsphere contact resonance as a single fit¬ 
ting parameter, and compare the fitted contact resonance 
with predictions based on microscale contact models. 

The metamaterial is composed of a monolayer of D = 
1.02 (im diameter polystyrene microspheres deposited on, 
and adhered to, the aluminum side of an aluminum- 
coated (100) silicon membrane of thickness 2 H = 1.3 pm, 
as shown in Fig. [T] The thickness of the membrane was 
measured with ellipsometry. The aluminum film is 50 nrn 
thick and serves as a medium to absorb the optical pump 
light. To fabricate the metamaterial, we utilize a self- 
assembly procedure, in which polystyrene microspheres 
are assembled at the air/water interface and then trans¬ 
ferred to the membrane, as shown in Fig.[l](a,b) [28!]. The 
resulting monolayer covers nearly the entire area of the 
membrane, which has dimensions of 4.8 mm x 4.8 mm. A 
representative microscope image of the monolayer pack¬ 
ing is shown in Fig. [He). 

The laser-induced transient grating technique [H, H3] 
used in this study is summarized as follows. The con¬ 
figuration of the beams and the sample is illustrated in 
Fig. He). Two optical pump beams derived from a pulsed 
laser source (532 nm wavelength, 430 ps pulse duration, 
and 1 kHz repetition rate) are overlapped at the alu¬ 
minum layer in the metamaterial, and form a periodic 
interference pattern. The pump spot has a 500 pm di¬ 
ameter at 1/e 2 intensity level. Absorption of the pump 
pulse light by the metamaterial induces a rapid ther¬ 
moelastic expansion, which induces counter-propagating 
acoustic waves with a wavelength equal to the optical in¬ 
terference pattern that is defined by the crossing angle 
of the beams. The acoustic wavelength is controlled by 
changing a phase mask pattern used to split the pump 
beam into +/ — 1 diffraction orders. The detection of 
the metamaterial acoustic response is accomplished with 



FIG. 2. Spectra of Fourier transform magnitudes where the 
black curve corresponds to the without-spheres case, the red 
curve to the with-spheres case, and each spectrum is normal¬ 
ized to its maximum amplitude. The markers denote the iden¬ 
tified peaks, which are plotted also in Fig. fusing the same 
markers and colors. The vertical lines denote the frequencies 
of the fitted microsphere contact resonance and spheroidal 
resonance. 


a quasi-cw probe beam (wavelength 514 nm and average 
power 10.7 mW at the sample) chopped to 50 ps pulses 
with an electro-optic modulator. The probe beam passes 
through the same set of optics as the pump beam and 
is focused at the center of the interference pattern to a 
spot of diameter 300 pm. Pump-induced surface ripples 
and refractive index variations caused by the propagat¬ 
ing acoustic waves (including contributions from the sil¬ 
icon membrane, aluminum film, and microspheres) lead 
to a time-dependent diffraction of the probe beam. The 
diffracted probe light is superimposed with an attenuated 
reference beam and is directed to a photodiode where it 
is recorded with an oscilloscope. 

Using this technique, we measured the acoustic re¬ 
sponse of the aluminum-coated membrane before the mi¬ 
crospheres were deposited (the “without-spheres” case). 
In this case, the pump pulses have an energy of 19 pJ 
and both pump and probe beams are incident on the 
aluminum-silicon interface in the metamaterial, such that 
they pass through the silicon membrane. As a result, 
the measured signal for the without-spheres case includes 
contributions from surface ripples from both the alu¬ 
minum film and the silicon membrane, and refractive 
index variations in the silicon. After the microsphere 
deposition, we measured the acoustic response of the 
metamaterial (the “with-spheres” case). In this case, 
both the pump and probe beams enter the metamate¬ 
rial from the opposite side, such that they are incident 
on the aluminum-microsphere interface, and the pump 
pulse energy is reduced to 5.6 pJ. Measuring the side 
with the microspher es p rovides enhanced sensitivity to 
microsphere motion [161] . in addition to signal contribu¬ 
tions from surface ripples of the aluminum film. 
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Figure [2] shows the normalized Fourier spectra of the 
acoustic oscillations corresponding to the without- and 
with-spheres cases at k = 0.78 /xm _1 . In the without- 
spheres case, two peaks are clearly observed. We identify 
the low frequency peak at 270 MHz as the Aq mode of 
the membrane, and the second peak at 1040 MHz, as the 
So mode. Because of the presence of the aluminum film, 
we predominantly excite antisymmetric modes, as can be 
seen by the relative amplitude of the Aq and So modes 
[29| . We measure the dispersion of these two modes by 
varying the acoustic wavelength, and plotting the peaks 
identified in the Fourier spectra, as shown in Fig. [3] The 
measured spectra and their identified peaks for all wave¬ 
lengths can be found in the Supplementary Information 
[301 ■ Using a silicon density of p m =2.33 g/cm 3 and typ¬ 
ical wave speeds in silicon [31] of cl = 9133 m/s (longi¬ 
tudinal) and Ct = 5844 m/s (transverse), corresponding 
to a propagation in the [110] direction with a transverse 
displacement along [001], we calculate the theoretical dis¬ 
persion curves for the Aq and Sq modes in the membrane 
using an isotropic model [32}. We find good agreement 
between our measurements and the calculated dispersion 
curves as illustrated in Fig. [3] We note that the mea¬ 
sured frequencies are slightly lower than predicted, which 
we attribute to our use of an isotropic model instead of 
the fully anisotropic model [33]. Calculations account¬ 
ing for the presence of the aluminum layer showed that 
the variation in the Aq and Sq dispersion curves due to 
the presence of the aluminum does not exceed 1% M- 
A clear difference between the with- and without-spheres 
cases can be seen by comparing the two spectra in Fig. [2] 
In the with-spheres case, each of the Aq and Sq modes 
identified in the without-spheres case are surrounded by 
two peaks. In this case, we also see an additional peak 
near 70 MHz. We observe this peak to be present in 
both the without- and with spheres cases for multiple 
wavelengths [30], and we do not consider it further for 
this study. 

The measured dispersion curves for the with-spheres 
case are plotted in Fig. [3] The with-spheres dispersion 
curves reveal “avoided crossing” [35] behavior between 
the Aq branch and the microsphere contact resonance. 
The lowest frequency branch of the with-spheres disper¬ 
sion curve follows the without-spheres Aq branch from 
the origin and then diverges to approach a horizontal 
asymptote near the contact resonance frequency. Above 
the resonance, a second branch follows the Aq branch at 
high wave vector magnitudes, but again diverges in the 
avoided crossing region to approach a horizontal asymp¬ 
tote just above the contact resonance. Measured peaks 
corresponding to dilatational motion closely match each 
other for the with- and without-spheres cases. The dis¬ 
persion curves also show a flat branch at f s = 1030 MHz 
that intersects with the Sq branch, which we attribute to 
microsphere spheroidal resonance. The spheroidal reso¬ 
nance frequency is obtained by averaging frequencies of 
the peaks identified along this branch. We note the pres¬ 
ence of two peaks at k = 0.78 /im -1 . This suggests either 


a hybridization between the between the Sq branch and 
the microsphere spheroidal resonance, or simultaneous 
detection of the Sq peak and the spheroidal resonance. 



FIG. 3. Dispersion relations. The red and black markers 
are the measured frequency peaks for the with- and without- 
spheres cases, respectively. The solid red line is the disper¬ 
sion calculated using our model. The dashed black line corre¬ 
sponds to the calculated So modes and the dashed blue line 
to Ao modes in the without-spheres case. The black dash- 
dot lines correspond to bulk waves in silicon. The horizontal 
lines denote the frequencies of the fitted microsphere contact 
resonance and identified spheroidal resonance. 

We compare the measured spheroidal resonance fre¬ 
quency with the frequency of the lowest frequency 
spheroidal mode /(1,2) of a free homogeneous isotropic 
polystyrene sphere, where the first index denotes the 
mode number and the second the harmonic number. Ap¬ 
plying the elastic equation for spheroidal modes [36|, we 
calculate a frequency of /(1, 2) = 1030 MHz using stan¬ 
dard longitudinal and shear velocities for poly styrene 
of Cl, p = 2350 m/s and ct, p = 1200 m/s [lCj, and a 
polystyrene density of p s = 1.06 g/cm 3 as provided by 
the microsphere manufacturer (Corpuscular, Inc.). As 
the calculated frequency is in close agreement with the 
measured frequency, we utilize velocities Cl, p and ct, p 
and density p s to solve for the polystyrene microspheres’ 
elastic modulus E s = 4.04 GPa and Poisson’s ratio v s 
= 0.32. The simultaneous identification of both the 
spheroidal resonance and the contact resonance offers an 
opportunity to study the variation of each within the 
same microsphere array, as has previously been demon¬ 
strated for isolated nanospheres [371 ]. 

To study Lamb wave propagation in our metamaterial, 
we adopt a similar approach as the one used previously 
to describe the interaction of a contact resonance of mi¬ 
crospheres with Rayleigh SAWs 0- We model the mi¬ 
crospheres as an array of linear surface oscillators that 
move in the Z direction, with mass m and linearized 
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normal contact stiffness K attached to the top surface 
of the thin silicon membrane, as illustrated in Fig. [ljd). 
The equation of motion for the surface oscillator can be 
written as mZ + K(Z — u z ^h) = 0, where u z h is the 
displacement of the substrate surface, Z is the displace¬ 
ment of the oscillator, and the microsphere mass m is 
calculated from the density p s . Because the frequency 
of the lowest frequency spheroidal mode is much higher 
than the contact resonance frequency and the elastic de¬ 
formation highly localized, we describe the entire mass 
of the microsphere moving in a rigid-body-likc motion at 
frequencies near the contact resonance. Since the acous¬ 
tic wavelength is much larger than the sphere size, we 
use an effective medium approach and approximate the 
average normal stress at the surface as the force exerted 
by the microsphere divided by the area of a unit cell [l6| . 
By applying boundary conditions corresponding to the 
average normal stress at one surface of the membrane in¬ 
stead of stress-free boundary conditions, we obtain the 
following dispersion relation for +x-propagating flexural 
and dilatational Lamb waves in the metamaterial [30]: 

2 4 

C ~~2 - 1 )DiD 2 = [ s inh (q L kH) sinh [q T kH)Di 

UJ q Zi P Crj-i 

+ cosh (q^kH) cosh ( qTkH)D 2 ], 

( 1 ) 

where qL = i/l — w 2 / (k 2 c\) , qx = \Jl — w 2 /(fc 2 c^), 
w c = 27r/ c = \JK/m is the angular frequency of the con¬ 
tact resonance, ui is the angular frequency of Lamb waves 
propagating through the metamaterial, A = y/3D 2 /2 is 
the area of a primitive unit cell in the hexagonally packed 
monolayer, and D 12 are determinants whose zeros com¬ 
prise the dispersion relations of uncoupled, purely flexu¬ 
ral and dilatational modes (respectively) in a membrane 
with stress-free boundary conditions |30j, |32|. The term 
in parentheses on the left-hand side of Eq. [T| describes 
the resonance of the oscillators. The right-hand side is 
responsible for the coupling between the oscillators, and 
flexural and dilatational Lamb waves. 

By taking the frequency of the contact resonance as 
the only fitting parameter, and using least squares mini¬ 
mization between the calculated and measured curves, 
we find the frequency of the contact resonance to be 
f c = 200 MHz. We plot the theoretical dispersion curves 
calculated using the fitted contact resonance frequency 
in Fig. [3] The theoretical dispersion curves show cou¬ 
pling between the flexural and dilatational modes due 
to the presence of the spheres adhered to one side of the 
membrane. In contrast to the coupling between Rayleigh 
SAWs and the microsphere contact resonance [l(j, the 
middle branch (as denoted in Fig. [1 does not stop at 
the line corresponding to transverse polarized bulk waves, 
and remains non-leaky due to the confinement of the 
membrane. This results in a second hybridization be¬ 
tween the middle and the Sq branch. The inset in Fig. [3] 
shows a closer view of this intersection, which is smaller 
than the avoided crossing near the Aq branch. This 


demonstrates stronger coupling between flexural motion 
and the contact resonance than with dilatational motion. 
This is a result of flexural modes having large out-of- 
plane displacement that couples with the vertical mo¬ 
tion of the surface oscillators, in contrast to dilatational 
modes with predominantly in-plane displacements. 

In Fig. [2] we observe that the largest signal contribu¬ 
tion comes from the peak nearest the contact resonance 
frequency. This is the case for each of the peaks along 
the flat part of the lower branch, which is consistent with 
previous observations [16] and is due to large signal con¬ 
tributions from the microsphere motion, as the oscillator 
equation of motion predicts large sphere oscillations and 
small surface displacements near the contact resonance. 
We observe a different behavior along the flat part of the 
middle branch [30l |. where the peaks closest to the mi- 
crosphere contact resonance are of small relative ampli¬ 
tude. This difference is also consistent with predictions 
from the oscillator equation of motion, as the flat part of 
the middle branch is further from the microsphere con¬ 
tact resonance as compared to the flat part of the lower 
branch. We also observe wider peak widths for peaks 
near the microsphere resonances despite low group veloc¬ 
ities. Possible causes include scattering due to disorder 
in the monolayer and inhomogeneous broadening due to 
variation in the contact stiffnesses. The latter includes 
both particle-substrate and particle-particle stiffnesses, 
although at the long wavelengths measured in this study, 
we do not expect significant interparticle effects. 

We estimate the frequency of the microsphere con¬ 
tact resonance using the Derjaguin-Muller-Toporov 
(DMT) adhesive contact model [la, [38] and compare 
this estimate with the fitted contact resonance fre¬ 
quency. We obtain a contact resonance / c ,dmt = 
(1/27T = 108 MHz, where Kdmt is the 
stiffness of the contact. To obtain the stiffness, we lin¬ 
earize the DMT contact model around its equilibrium 
position [30i ], as no anharmonic behaviors are observed 
for these pump powers, and we estimate the displace¬ 
ment of the surface to be much smaller than the equi¬ 
librium overlap of the microspheres [3Cj. As for previous 
measurements of 1.08 /im diameter silica microspheres 
on a thick aluminum coated fused silica substrate [l6[ . 
we find that the measured contact frequency is much 
larger than the frequency predicted utilizing the DMT 
contact model. The use of alternative contact models re¬ 
sults in little variation in the predicted contact resonance 
frequency [30]. Instead, we suggest that other uncertain¬ 
ties in the modeling of the contact may be the cause for 
the discrepancy between the estimated and the measured 
values of the contact resonance frequency. 

We have presented the realization of a self-assembled, 
locally resonant metamaterial for Lamb waves composed 
of a monolayer of microspheres adhered to a thin elas¬ 
tic membrane. With their high characteristic frequen¬ 
cies and small length scales, these metamaterials hold 
promise for the development of new types of Lamb wave- 
based devices. Since the response is sensitive to the state 
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of the contact, these metamaterials may have potential 
applications as sensors for humidity, temperature, and 
micro/nanoscale material properties, and serve as a plat¬ 
form for the exploration of microscale contact mechanics. 
Because of the scalability enabled by self-assembly, this 
type of metamaterial may enable future studies that ex¬ 
plore the interaction of local mechanical resonances with 
even higher-frequency phonon transport. Finally, the use 
of self-assembly approaches may also offer significant po¬ 


tential advantages for producing acoustic metamaterials 
inexpensively and on a large scale. 
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